A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
02(12.03.1; 12.03.3; 03.13.6) 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



Peak statistics on COBE maps 

R. Fabbri^ and S. Torres^ 

^ Dipartimento di Fisica deirUniversita, Sezione di Fisica Superiore, Via S. Marta 3, 1-50139 Firenze, Italy 
Observatorio Astronomico Nacional, Universidad Nacional de Colombia, Bogota, Colombia 

Received ; accepted 



Abstract. We perform the stastistics of temperature 
maxima and minima in COBE-DMR 2-year maps. For 
power-law spectra the surface distribution of peaks implies 
an amplitude consistent with more conventional analyses 
of COBE data (for instance, we get Qrms-ps = 17± 3 /iK 
for a spectral index n = 1), but not with the measured 
quadrupole Qrms = 6 ± 3 /iK. This provides further sup- 
port for the existence an infrared cutoff in the cosmic spec- 
trum. 
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1. Introduction 

The COBE experiment (Bennet et al. 1992, 1994; Smoot 
et al. 1992) has stimulated a considerable amount of work 
on cosmic structure. Current tests usually exploit the an- 
gular correlation function and several harmonic ampli- 
tudes of the sky temperature field (see e.g, Adams et al. 
1992; Kashhnsky 1992; Efstathiou et al. 1992; Kofmann 
et al. 1993; Gorski et al. 1994). However, several more 
tests have been suggested over the years for a thorough 
investigation of the properties of the anisotropy field of 
the cosmic background radiation. These often involve the 
distribution and features of hot and cold spots, which can 
provide useful checks of the Gaussian nature of the fluc- 
tuations (Sazhin 1985; Bond & Efstathiou 1987; Coles & 
Barrow 1987; Coles 1988; Martinez-Gonzalez & Sanz 1989; 
Gott et al. 1990). Measurable quantities include the num- 
ber of spots Aiso defined by isotemperature contours, the 
spot boundary curvature or genus G, the spot excursion 
area, and so on. A first analysis of COBE-DMR maps 
along these lines has been performed by Torres (1994). 
A definition of spots independent of isotemperature con- 
tours considers local maxima and minima of temperature 
(Bond & Efstathiou 1987, Vittorio & Juszkiewicz 1987), 



and is thereby not connected to the topological features of 
spots. Differing from Aiso, the dependence of the number 
of (positive and/or negative) peaks Apeak on threshold is 
not universal for Gaussian fields. This latter approach was 
adopted by Fabbri & Natale (1993, 1995) in studies of the 
2-dimensional distribution of extragalactic IRAS sources, 
but has not yet been applied to the cosmic background 
radiation. 

In this work we analyze the statistics of local maxima 
and minima in COBE-DMR 2-year maps. We found that 
in this kind of analysis the detector noise must be taken 
into account very carefully since Apeak is sensitive also to 
high order harmonics where noise dominates (cf. Fabbri 
1992). However, due to a highly nonlinear dependence of 
Apeak on the harmonic strengths, the presence of struc- 
tured signals in COBE maps reduces its value below the 
level measured in pure noise maps. (The identification of 
genuine peaks in the radiation temperature is not required 
at all in our analysis.) We find that the distributions of 
positive and negative peaks are mutually consistent, and 
the results from this statistics agree with those of earlier 
tests. Therefore, we find no evidence for non- Gaussian fea- 
tures in the fluctuations. More precisely, fitting Gaussian 
power-law models of cosmic structure to the peak distri- 
bution we recover a clear anticorrelation between the spec- 
tral index n and the predicted rms quadrupole Qrms-PS 
(Seljak & Bertschinger 1993; Smoot et al. 1994): We get 
Qims-ps = 17 ± 3 /iK for n = 1 and 14 ± 3 /iK for 
n = 1.5, where the error bars include uncertainties deriv- 
ing from the treament of noise as well as cosmic variance. 
These numbers, altough they agree with previous eval- 
uations of the quadrupole from higher order harmonics, 
are not consistent with its direct determination providing 
Qrms = 6±3 /iK (Bennet et al. 1994). So the recently dis- 
covered discrepancy is confirmed by the properties of the 
peak distribution, which depend on the harmonic content 
of the angular distribution up to ~ 50. 
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2. The peak number test 

We analyzed the 2-year 53(A+B) DMR maps processed 
with a 2.9° smoothing (Wright ct al. 1993) and dipolc sub- 
traction. Considering the Northern and Southern hemi- 
spheres separately, we constructed two pole-centered 
maps, each containing 12892 pixels, using the coordinate 
transformation 9i = 2sin[i (| — |6|)] , = Z, with b 
and I the Galactic coordinates. After masking low Galac- 
tic latitudes, |6| < 20°, wc were left with 8412 pixels per 
map. We then looked for temperature peaks using the al- 
gorithms of Fabbri & Natale (1993, 1995). 

Table 1 gives the no-threshold numbers of peaks, both 
actually detected and extrapolated to the entire sky and 
to the North and South hemispheres (2nd and 3rd column, 
respectively). Figure 1 reports the extrapolated numbers 
vs. a threshold factor ly. This is the peak height normal- 
ized to the sky rms fluctuation (75(0); for cold spots, 
^poak gives the number of minima below — C2(0)z/. For 
distributions of only positive or negative peaks the sta- 
tistical errors at 1-sigma confidence levels are evaluated 
as (A''pcak/./u) ^ ! with fu = 0.652 the unmasked fraction 
of the sky. Within such error limits, we find no signifi- 
cant difference between the distributions of positive and 
negative peaks; this result provides support for the Gaus- 
sian nature of cosmic perturbations. In the figure we also 
report the full-sky average number of positive and nega- 
tive peaks. This average will be compared with theoretical 
models below, because of the smaller (by a factor v^) rel- 
ative error. 

For the interpretation of the above data, we must con- 
sider that both cosmic signal and noise contribute to the 
temperature field T{9,(j)) = J2^^fmYem{d,4>)- The expec- 
tation values of the numbers of local maxima and minima 
are determined by the variances a| = J2m (I 1^)- ^'^^ 
a Gaussian field the full-sky A^peak is given by 
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where 

z,i = (3-27*')-*, U2=[2(l-7*')]"*, 



(1) 



(2) 



and the properties of the anisotropy field are summarized 
by the parameters Q* and 7* (Bond & Efstathiou 1987; 
Fabbri 1992) given by 



Fig. 1. The peak number A'^peak vs. threshold in COBE-DMR 
maps. Open squares (triangles) denote positive (negative) 

peaks. The four lowest curves give the peak numbers ex- 
trapolated to the Northern (dotted Unes) and Southern 
(dash-dotted) halves of the sky. The remaining curves refer 
to the full sky. Filled circles give the average numbers of pos- 
itive and negative peaks. The full line describes a power-law 
model with n = 1 and Qrms-ps = 18.5 /^K, corresponding to 
the fit procedure S31 in Table 2 
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and 
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Y,{e-l)e{£+l){£ + 2)W!aj 
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Here We are form factors taking into account beam shape 
and any additional smearing effect. 

Assuming that signal and noise are uncorrelatcd, their 
contributions to aj add up in quadrature in the (A-|-B) 
maps, a| = -|- a^^, and we need to determine a^^ in 
an independent way. This can be achieved by means of 
the (A— B) maps, which can however be used to directly 
derive the coefficients = W^a^f rather than the a^^. 
Figure 2 reports the results of a harmonic best fit up to 
£ = 30 executed on the entire celestial sphere. Increasing 
the number of harmonics, best fits (as well as direct in- 
tegration by means of a^,„ = J T (9 , (ti)Yf*j^{9 , (j))dil) tend 
to overrate large-i? amplitudes. We checked for this effect 
by considering the peak statistics for noise maps. Figure 
3 gives the numbers of maxima and minima detected in 
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Table 1. No-threshold peak numbers 



Peak set 


{A + B) Maps 
|fe| > 20° 


{A + B) Maps 
Extrapolated^ 


{A - B) Maps 
No mask 


Hot, North 


50 


76 6 ± 10 8 


95 


Cold, North 


47 


72.0 ± 10.5 


99 


Hot, South 


48 


73.6 ± 10.6 


105 


Cold, South 


53 


81.2 ± 11.2 


97 


Hot, total 


98 


150.2 ± 15.2 


200 


Cold, total 


100 


153.3 ±15.3 


196 


Hot /Cold Average 


99 


151.7 ± 10.8 


198 



To full sky or hemispheres, including masked regions 



Fig. 2. The noise harmonic coefRcients A^g (containing the 
form factors Wf) and the beam shape coefficients Gi vs. (.. 
Filled circles give as obtained from the noise map, and 
the dotted line represents the best fit curve corresponding to 
o-N = 0.936° and Cn = 5.74 x 10"^ mK^. Crosses give the 
experimental Gf, compared to Gaussian approximations with 
dispersion 2.9° (full line) and 3.1° (dash-dotted) 

the full-sky noise maps vs. threshold v. If using Eq.s (1-4) 
we generate the peak statistics from a set of 50 harmonics 
with best-fitted amplitudes, we get a low-threshold excess 
of about 30 peaks with respect to the data in Fig. 2; the 
discrepancy increases with the number of harmonics. Since 
accurate theoretical calculations of peak statistics require 
at least ^^50 harmonics, we tried to fit an analytic form to 
a more limited set of A^^. A satisfactory choice is 




Fig. 3. The peak number iVpoak vs. threshold in COBE-DMR 
noise maps. Open squares (triangles) denote positive (nega- 
tive) peaks. The four lowest curves give the peak numbers 
extrapolated to the Northern (dotted lines) and Southern 
(dash-dotted) halves of the sky. The remaining curves refer 
to the full sky. Filled circles give the average statistics of posi- 
tive and negative peaks. The full line represents the best fitting 
function 

However the large error bars on individual make a 2- 
parameter fit for the function (^ not very useful. A more 
convenient solution to this problem is to fit the parameter 
ctn directly on the noise-generated peak distribution; us- 
ing the average distribution of maxima and minima (rep- 
resented by the filled circles in Figure 3) and 99 harmonics 
the result is ctn = 0.936° ± 0.071°. (We checked however 
that a set of ~ 50-^60 harmonics would provide a sufficient 
accuracy.) The fit turns out to be even "too" good, provid- 
ing Xmin = 1-3 against 9 degrees of freedom: This means 
that data points at different thresholds are not uncorre- 
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Fig. 4. The acceptance regions in the (n, Qrms-ps) plane at 
1- and 2-sigma confidence levels. Contours are computed with 
procedures S31 (full lines) and F31 (dotted). Both procedures 
= 3.1°, (JN = 0.936°, and Cn = 5.74 • 10~^ mK^ 



assume as 

The former uses the measured A^g up to 



30. 



lated. The value found for ctn is quite independent of Cn, 
since peak statistics does not depend on the overall ampli- 
tude of anisotropics as shown by Eq.s (3, 4) .We then fitted 
the amplitude parameter of Eq. (^) on the reported set of 
A^^ with ^ = 2 H- 30, getting Cn = (5.74 ± 0.66) x IQ-^ 
mK^. The distribution described bu the full line in Fig. 3 
was generated using the function (^) in the entire range 
i = 2 99. We checked that the peak distribution of the 
(A-B) maps is also recovered to a good accuracy building 
up the harmonic spectrum with the "measured" A^^^ up 
to = 30 and the function (||) at higher £. 

We tested power-law spectra against the peak statis- 
tics in Fig. 1 with the following procedure. For each spec- 
trum, labelled by the spectral index n and the predicted 
quadrupole Qrms-PS = <iS2/V^, we generated the theo- 
retical expectation values Og^ and then 



W, 



2„2 



exp 



2 



A 



(6) 



Here Gi denote the measured beam-shape coefficients of 
COBE-DMR. These have reported by Wright et al. (1994) 
up to ^ = 50. A more extensive set of Gi up to £ = 99 has 
been provided to us by Kogut (private communication). 
The exponential factor in Eq. (^ takes into account the 
2.9° smearing on the map as well as the additional smear- 
ing due to orbital motion, so that our best estimate is 
as = 3.1°. Note that as should not be confused with 
the approximate Gaussian beamwidth which was used in 
many computations, but not in the present work. Also, 
it is quite different from the phenomenological parameter 
aN of Eq. (H). The peak distributions calculated from Eq.s 
(1-4) are tested against the averaged distribution of posi- 



tive and negative peaks. We should notice that such the- 
oretical distribution are affected by several sources of un- 
certainty. Cosmic variance (White et al. 1993) affects the 
cosmic-structure contribution to the harmonic coefficients 
aim, and a quite similar variance pertains to noise; in fact, 
these two effects are described by identical equations as- 
suming that both are Gaussian processes. Performing a 
limited number of simulations for superpositions of Gaus- 
sian signal and noise (with fixed aj , the expectation values 
of harmonic strengths) , we found that the probability dis- 
tribution of iVpcak at a given threshold is roughly consis- 

tent with a Poisson distribution having a width (iVpcak) ^ • 
However, in our case the Af^^ themselves are not fixed. 
When we use Eq. (||) their errors at different £ are cor- 
related due to the uncertainties on Cn and ctn; it can be 
shown that the corresponding contribution to the uncer- 
tainty on the predicted A'pcak is given by 



var(iVpcak) 




d (A^poak) \ var(CN) 



dA^g 



C2 



E 



A 



d (A^peak) 
9^N« 



(TNvar((TN)- 



(7) 



When we use the measured A^^, a more familiar equation 



var(Afpcak) = 



applies. 



(a(A^peak) /a^y var(A^,) 
For each computed model we calculated the corresponding 
combining quadratically the above errors with with the 
experimental error bars reported in Fig.l. This procedure 
allowed us to avoid a more extensive use of simulations. 

Figure 4 gives the allowed regions at 1 and 2 sigmas in 
the (n, Qrms-ps) plane. Full lines give the contours calcu- 
lated for as = 3.1°, using the of Fig. 2 up to £ = 30, 
and the analytic form (|^) with the optimal values of ctn 
and Cn at higher £. (This procedure is referred to as S31 
in the Figure and in Table 2.) Contours calculated for 
as — 2.9° (case S29, not reported in the Figure) would 
be hardly distinguishable from them. Using the analytic 
form (H) in the entire range £ = 1 99 (case F31) the con- 
tours are slightly displaced to smaller values of Qrms-ps 
(by - 2 for n = 1). 

The parameters n and Qrms-ps are clearly anticorre- 
lated, as already found from analyses of harmonic am- 
plitudes and the angular correlation function (Seljak & 
Bertschinger 1993; Torres et al. 1994; Smoot et a. 1994). 
Minimizing for fixed n we identify a straight line in the 
(n, Qrms-ps) plane. 



Or 



-PS = a + hn, 



(8) 



with the coefficient values listed in Table 2. From the 
above results, taking into account differences arising 
from the S31 and F31 procedures, we can conclude that 
Orms-ps = 17 ± 3 for n = 1, and Qmis-ps = 14 ± 3 
/iK for n — 1.5. These numbers agree very well with the 
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Table 2. Coefficients of tlie n — Qrms-ps regression defined by 
Eq. § 



Procedure 






S3lt 


24.71 ±0.65 


-6.23 ±0.45 


S29 t 


24.54 ± 0.66 


-6.16 ±0.46 


F31* 


21.47 ±0.60 


-5.45 ±0.42 



'ctn = 3.1°, using noise tiarmonic amplitudes from Fig. (2) for 
£ < 30. 
*crN = 2.9° 

* Noise described by Eq. (0) for any I 



most likely quadrupole Qrms-PS derived from higher order 
multipoles, but not with the quadrupole rms fluctuation 
of 6±3 /iK directly fitted on two-year data; see Bennet et 
al. (1994) for a discussion of this discrepancy, and Jing 
& Fang (1994) for a possible explanation in terms of an 
infrared cutoff in the spectrum. 
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